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Abstract — Reasoned by its dynamical behavior, the memristor 
enables a lot of new applications in analog circuit design. Since 
some realizations are shown (e.g. 2007 by Hewlett Packard), the 
development of applications with memristors becomes more and 
more interesting. Whereas most of the research was done in 
the direction of memristor applications in neural networks and 
storage devices, less publications deal with practical applications 
of analog memristive circuits. But this topic is also promising 
further applications. Therefore, this article proposes a frequency 
dependent rectifier memristor bridge for different purposes (e.g. 
using as a programmable synaptic membrane voltage generator 
for Spike-Time-Dependent-Plasticity) and describes the circuit 
theory. In this context it is shown that the Picard Iteration is 
one possibility to solve the system of nonlinear state equations 
of memristor circuits analytically. An intuitive picture of how a 
memristor works in a network in general is given as well. In 
this context some research on the dynamical behavior of a HP 
memristor should be done. 



I. Introduction 

THE usage of memristors in analog circuit design enables 
new applications. In (TJ, an ADC consisting of memris- 
tors has been proposed. Another applications are an automatic 
gain control circuit [2| .programmable analog circuits Q, 
an electrical potentiometer Q or oscillators EJ, 0. In |6) 
memristors are used for basic arithmetic operations. 

In this paper a frequency dependent rectifier memristor 
bridge is presented. Therefore, a general description of mem- 
ristive systems will be given first. Then the HP memristor 
is presented and some analysis about its dynamical behavior 



will be shown in III and in IV properties of memristors are 



given. In[VJ the memristor bridge is presented, which is solved 
analytically by Picard Iteration ir VI Finally, the results of the 



present investigations are summarized. 

In General, a memristive system is described by the two 
equations Q 

V(t)=M(x,I,t)-I(t) (1) 



dx 



(2) 



where V(t) is the applied voltage and I(t) is the current 
through the device. M(x,I,t) is called memristance and relates 
voltage and current. Thus, similar to a linear resistor its 
dimension is Ohm. x is a vector of the dimension n which 
consists of the internal state variables. The dynamics of this 
vector is specified by the state equation [2] 



II. HP MEMRISTOR 

The HP memristor was realized in 2007 by a team of 
researchers of Hewlett Packard around Stanley Williams [8|, 
|9l . While this memristor is passive, latest realization of a 
memristor is an active one on the base of niobium oxide [10|. 
The HP memristor is made of a titanium dioxide layer which 
is located between two platinum electrodes. This layer is in the 
dimension of several nanometers and if there is an oxygen dis- 
bonding, its conductance will rise instantaneously. However, 
without doping, the layer behaves as an isolator. The area of 
oxygen dis-bonding will be referred to space-charge region and 
changes its extension if an electrical field is applied. This is 
done by drifting the charge carriers. The smaller the insulating 
layer the higher the conductance of the memristor. Also, the 
tunnel effect plays a crucial role. Without an external influence 
there is no change in the extension of the space-charge region. 
The internal state x is the normed extent of the space-charge 
region and can be described by the equation 



w 



x := — , <x < 1, x £ ' 



(3) 



where CO is the absolute extent of the space-charge region and 
D is the absolute extent of the titanium dioxide layer. The 
memristance can be described by the following equation |9|: 



M{x) = R 



R ff(l-x) 



(4) 



R on is the resistance of the maximum conducting state and 
R Q ff represents the opposite case. The vector of internal state 
of the HP-memristor is one dimensional. That is why the scalar 
notation is used. The state equation is 



dx 
df 



■■f(x,I,t) 



D 2 



T M (0 



(5) 



where jj, v is the oxygen vacancy mobility and is the 

current through the device. The following quantizations are 
used for memristors in this paper. 



Ron : 

Roff: 



100 Q., 
I6k£l, 



D 



10" 
10" 



-14 



s-V 



These values are given as an example by the Hewlett Packard 
Team. The ratio of maximum and minimum value of the 
memristance can be as large as e.g. 160. 
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A. Window function 

The extension of space-charge region is physically limited. 
Therefore, a window function shall be established to take ac- 
count of saturation which is needed for a more realistic model. 
Several shapes of window functions are possible e.g. ifTTl . In 
the following, the window function 



if < x < 1 
else 



(6) 



is used fl2l . The current mathematical description of this 
window function shows a weakness in numerical calculations. 
Excluding external influences, in the case of {x < or x > 1 } 
the window function will be 0. Therefore, a state change of the 
system will not be possible. But considering Vm, the voltage 
over the memristor, for the two cases x < and Vm > or 
x > 1 and Vm < 0, a state change should be possible. The 
window function will result in: 



f w (x,V M ) = 



(4-x- (1 -*))*, ife<x<l-e 

or x < e and Vm > 
or x > 1 — e and V M < 
0, else 

(7) 

A better stability for numerical simulations is achieved by in- 
troducing e and therefore an extended restriction of definition 
area. For all simulations in this paper e equals 0.02. 

B. Polarity of the HP memristor 

An agreement concerning the memristor polarity will be 
introduced. The extension of the space-charge region is going 
together with a decrease of the memristance and vice versa. 
The polarity determines the behavior of the memristor. Ac- 
cording to the direction of the applied field, if the memristance 
decreases the memristor will be "forward biased". The other 
direction should be called "reverse biased". Because of the 
minimum resistance and the time lag, both definitions are main 
intuitive. For a memristor the algebraic sign of Vm equals that 
of Im and the polarity just depends on it. Including the window 
function and the algebraic signs, Eq. [5] extends to 

dx _ ji v ■ R 
df ~ D 2 

if the memristor is forward (reverse) biased. 

III. Time response of a HP memristor 

A. Change of the internal state 

The calculations of this part will be done by neglecting the 
window function. The aim is to investigate the time behavior 
of an AC sine current source Is which is connected in series 
with a HP memristor. Regarding the state equation |5]l 

dx _ fi v ■ R 
df D 2 _ 

Is=Im 



■fw(*, + (-)V M )-( + (-)I M (0) (8) 



•/o'sin(co-f) 



(9) 



dividing of the variables on both sides 



pxq+Ax 



dx = const ■ Iq 



t +At 



sin((0-t')dt', Af >0 (10) 
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(a) current source: Iq = 0.0005A, (b) current source: /o = 0.0022A, 
voltage source: Vq = IV. voltage source: Vq = 30V. 

Fig. 1. Dashed line: supplied sine current source, solid line: supplied sine 
voltage source, f = 30Hz, (a) Maximal state change Ax maI in dependency of 
xq, At = E, (b) M(x) in dependency of t for two different initial states and two 
source types. 



and performing integration leads to 

Ax = - — ™ HSt ■ ^ cos(co ■ (to + Af)) - cos(co - to)) . (11) 

whereas xo :=x(to). At (to + Af) the extension of the space- 
charge region is (xo + Ax). In the case of a memristor which 
is connected in series with a sine current source, the change 
of the space charge region Ax does not depend on the initial 
state xo (see Eq. [TTJ. There will be an other result if there is 
a sine voltage source instead. In this case the state equation is 



dx 
df 



Mv • R o 
D 2 



Vb-sin(ow) 
M(x) 



(12) 



and because of M(x) it leads by integration to a quadratic 
equation. The relevant solution of this quadratic equation for 
Ax with to = is 



Ax=(r-x )-J(r- x ) +b ■ (cos(a) • At) - 1) (13) 



whereas r : 



and b : 



2-VQ-const 



The second solu- 



R off- R on ' ra-(R ff- R oL, 

tion is irrelevant because for At = the change of the space 
charge region Ax has to be zero. Reasoned by the disregard of 
the window function, the usage of equation ( p"3] > is bounded 
because the solution has to be within the physical limits. 
But the bottom line is that Ax depends on the initial state 
xo, if there is a voltage source in series with the memristor. 
As shown before, for a supplied current source it does not. 
Fig. [TJa) illustrates these insight. Since the memristance is 
determined by the internal state, the change of the memristance 
behaves similarly. As you can see in figure [TJb) for a sine 
voltage source the change of the memristance depends strongly 
on the internal memristance, however for a sine current source 
it does not. Note, in the case of a sinusoidal signal, the 
maximum change of the space charge region Ax muv will be 
reached at the end of one half period. This is because the 
direction of the change of the space charge region depends on 
the algebraic sign of the external signal. 
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B. Highest frequency for fully state change 

To estimate the dynamical behavior of the HP-memristor in 
circuits the frequency f cut for sinusoidal signals should be in- 
troduced. This is the highest frequency for which a memristor 
will be able to change from lowest to highest memristance or 
vice versa. It means that this state change completes exactly 
at the end of one half period (At cut = 0.5 • T cut ). Note, for 
the used memristor model and same conditions, the required 
time for changing from a state A to a state B is the same 
as for the reverse process. f cut depends on the amplitude 
of the supplied source. Neglecting the window function the 
calculations should be performed for a single memristor in 
series with a supplied sine voltage source. This example allows 
rough estimates for circuits which are more complex. Taking 
into account that CO equals 2-n- f and T equals 4. Integration 
of Eq. [13] and solving for / with t = leads to 



1 

% 



Vq ■ const 



+x • Ax) • (Ron - R off ) + Ax • R off 



(14) 



whereat f cut is direct proportional to the amplitude Vo. 
example: Considering the restricted domain of x the maxi- 
mal frequency for switching from highest memristance x(t = 
0) = 0.02 to lowest memristance x(t — -|) = 0.98 should be 
calculated. Therefore Ax is equal 0.96. Calculating by using 
Vq = 30V leads to f cu t ~ 12.35675//z. Testing by solving the 
state equation numerically leads to following results: 
x{t = ±j)= 0.98 for / = 12.35675 Hz 
x ( t= Tf)= 0-9424 for / = 12.4#z 

If the frequency is higher than f cut , saturation wouldn't be 
reached. Considering the window function, the result of solv- 
ing the state equation by simulation is x(t = ^j) = 0.98 
for /= 10.75 Hz. Therefore the calculations without window 
function are reasonable for rough estimates. Note, for / > f cut 
if the frequency is increasing the ratio of maximum and 
minimum value of the memristance will decrease. Note, for a 
supplied sine current source the calculation for f cut is also 
possible (Conversion of Eq. 



amplitude Iq 



30V 
R off 
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As an example, 
= 1.9mA, f cut is about 6.2Hz- 

IV. Properties 



for an 



To summarize, the value of memristance depends on the 
current load of the past 1131 . If there is no current flow, the 
internal state will be retained. From this it follows that a 
memristor acts like a nonvolatile memory, whereas the range 
of values is continuous [14], ifTBll . That is an interesting fact 
comparable to transistor memory technology. The difference 
between highest and lowest possible memristance is relative 
large [9]. Indeed in the true sense, the memristor is no switch 
but it could be used for switching operations. Therefore, 
depending on the direction and the benchmark the memristor 
passes higher potentials and locks for lower ones. Subject 
to the time shift the behavior of the memristor is frequency 
dependent. For lim/^oo it behaves like a linear resistor Q 
because the change of the internal state couldn't follow the 
rapid voltage change. For sufficiently small frequencies the 
nonlinearities are dominating whereas the time shift is direct 
proportional to the amplitude of the signal. 




V M i 




Vio 


(15) 


Vm2 


= VlQ 




(16) 


Vm3 


= v s - 


v 20 


(17) 


Vm4 


= V 2 Q 




(18) 


Vrl 


= v lQ - 


Via 


(19) 



Fig. 2. Schematic of the memristor bridge circuit. 



V. Memristor bridge circuit 

This circuit (Fig. [2]) contains four HP memristors, one AC 
voltage source and one load resistor. Per definition for a 
positive voltage, memristors Mi 4 are forward biased, while 
M23 are reverse biased. In [16] this circuit is also been 
presented, but there is a difference in application. While 
this paper deals with periodic signals and their specifics at 
different frequencies, in |fl6l pulses are used for synaptic 
weight programming. 



A. Mathematical describtion 

At this point the new notation for the memristance 

M(x n ) :=M„, n e N 

should be established. The voltages 

Num + M 2 ■ M3 • R L 



V 



L0 



Vs 



V20 :=V2I0(XI,X2,X3,X4) = V s • 



Den 

Num + Mi • M4 • Rl 
Den 



(20) 



(21) 



(22) 



are defined from the nodes to ground, whereas the notations 

Num : = (Mi + M 3 + R L ) • M 2 ■ M 4 (23) 
Den := (Mi + M 2 ) ■ (M 3 -M 4 + M 3 • R L + M 4 • R L ) 

+ Mi ■M 2 -(M 3 +M 4 ) . (24) 

are used for a better overview. The structure of the circuit 
implies a nonlinear system of differential equations of the 
fourth order. The four state equations are 

dx„ 



' Ron „ , , n , Vm„ 

hln 



n <4 



(25) 



For xi )4 the sign is "+" and for x 2 ^ the sign is negative. 

B. Simulation results and functionality 

Using these state equations, numerical simulations are pos- 
sible. From this point Vs should be a sine AC voltage source 
and R L = 1 k£l. 

Regarding Fig. [3] dependent on the frequency a qualitative 
difference for the voltage over the load Vrl is detectable. 
For low frequencies (represented by / = 0.5 Hz in the case 
of Vq — 30V) the output signal is almost exclusive positive, 
comparable with a rectifier circuit. Therefore, the frequency of 
the output signal is twice as high as the one of the input signal. 
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1 2 3 4 0.01 0.02 0.03 0.04 0.05 0.06 

Time (s) Time (s) 

(a) (b) 

Fig. 4. Memristances M\ and M3 over time t, sine AC voltage source, Vq = 
30V, x„(t = 0) = 0.5, (a) f = 0.5Hz, (b) f = 30 Hz. 



Reasoned by the time shift of the change of state, which is 
subject to a HP memristor, there are temporary short negative 
peaks. Using a higher excitation frequency (represented by / = 
30 Hz for Vo = 30 V) there is no rectifier function detectable. 

By means of figure [4] this frequency selective behavior can 
be better understood. For low frequencies (/ < f cut ), a com- 
plete state change will happen, the memristors are going into 
saturation. Beginning with an initial value M\ ,4 is decreasing, 
while M23 is increasing until saturation is reached. When the 
sine voltage source reaches its negative half period the process 
will be reversed until saturation is reached again and so on. 
During one period there is a change for the relational operator 
between both memristances. So every half period, the potential 
in node one is going to be higher than the potential in node 
two, which implies that Vrl is always going to be positive. 
For / > f cut , the memristors do not completely switch from 
lowest to highest memristance. It follows that they are not 
going into saturation and the memristors will just reach the 
initial state after one period. So there is no change for the 
relational operator between M\ .4 and M23 during one period 
and the relation of the memristances just depend on the initial 
conditions (See figure Qb)). For example, if the initial states 
are equal for all four memristors, then M\ 4 would be less than 
M23. So Vrl would be positive for the first half period too. But 
at the second half period, the Memristances M2 and M3 keep 
higher than M4 and Mi . So the potential at node "2" is higher 
than at node "1", which implies that Vrl is negative. There are 
two reasons why the absolute value for the amplitude of Vrl 




(a) lim/_>o (b) lim/ 



Fig. 5. Approximate circuit equivalents: (a) Graetz circuit for very low 
frequencies, (b) Wheatstone circuit for very high frequencies. 




(a) (b) 



Fig. 6. Vrl if dependency oft, supplied periodic square wave voltage source, 
V = 10 V, x„ [t = 0) = 0.5, (a) f = 0.5Hz, (b) f = 30H Z . 

decreases for increasing frequency. On the one hand for equal 
initial conditions the maximum potential difference between 
node "1" and "2" decreases by increasing the frequency. On 
the other hand, dependent on the initial conditions, there is 
may a higher voltage drops across the memristors. Note, the 
maximum of the amplitude of Vrl is not at t = ? . At this time 
the amplitude of the sine source begins to decrease, but the 
change of the memristor states is continuing, which has firstly 
more impact. 

The lower the frequencies the more the circuit behaves like 
a Graetz circuit. For very high frequencies the circuit is going 
to behave like a Wheatstone bridge. This can be proven by 
the voltage divider. For liiriR^oo it simply is 

Vrl = V10 - V 2Q = V S ■ ( - Mz - - M \ 4 ) , (26) 
\M\ +M2 M3 +M4 / 

and had to be zero for equal resistances, if the conclusion 
is true. Simulations illustrate that for e.g. / = 50kHz and 
x n (0) — 0.5 the amplitude of Vrl equals 0.5 mV, while the 
ratio between highest and lowest memristance of e.g. M{x{) 
equals 1.0002. 

Using a supply periodic square wave voltage source with 

V s =^«(sin(o)-f)). (27) 

is another interesting example. The results for this are shown 
in Fig. [6] A frequency dependent behavior is also detectable. 
For higher frequencies the voltage over the load resistor is 
serrated. For low frequencies this voltage is almost constant 
with negative peaks. The condition that the initial states of 
all memristors are equal leads to the in Fig. [6] presented 
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Fig. 7. (a) System supplied by a periodic square wave voltage source with 
V Q = 10 V, f = 30Hz, -vi (0) =*t(0) = 0.5 +e, .v 2 (0) = x 3 (0) = 0.5 - e. (a) 
Mi (x) and M3 (x) in dependency of t for different initial conditions, (b) Vrl in 
dependency of t for different initial conditions. 



behavior. But what happens if the circuit is supplied by a 
periodic square wave voltage and the initial states are not 
equal? As it is described in chapter |Hl| and shown in Fig. |7ja), 
for a supplied periodic voltage and higher frequencies the 
variation of memristance depends on the initial conditions. 
So the voltage curve changes from the saw tooth form. As 
represented in Fig. |7|b) it is possible to create a curve which 
is similar to synaptic membrane voltages which are used for 
Spike-Time-Dependent-Plasticity [17|. The amplitude of the 
voltage varies for different initial conditions, which could be 
set by a low frequency or DC voltage. 



C. Possible applications 

Depending on the frequency, there are two significant 
species. Fusing a common Graetz circuit and a Wheatstone 
circuit in one device is one application possibility. Replacing 
the independent source with a controlled one leads to a further 
application. The source is controlled by the load voltage and 
its frequency is variable. Keeping in mind, the amplitude of the 
output signal decreases by increasing the frequency. Because 
of that, a frequency controlled regulator is conceivable. If 
the amplitude exceeds a predefined value then the frequency 
should be increased to prevent a further amplification. Using 
the circuit as a saw tooth generator is also possible. Like 
mentioned before, the presented circuit also could be used 
as a programmable synaptic membrane voltage generator for 
Spike-Time-Dependent-Plasticity. 

VI. Picard Iteration 

At this point, the Picard Iteration should be introduced as a 
possibility to solve common memristive systems analytically. 
One advantage related to e.g. the Volterra-series expansion 
is that the Picard Iteration converges more rapidly. In this 
chapter, the application possibility of this Iteration should be 
shown on the memristor bridge circuit. In favor, the general 
description of the Picard Iteration should be given first: A 
nonlinear, dynamical system with 



dt 



f(t,x(t)), x(to) 



- xq 



(28) 



is given, whereas x(t) is Lipschitz continuous and an element 
of a Banach space. Then, the iteration is given by 



Xq 



+ f f(t ,x {k \t'))-dt , fe{f ,fo + e} 

Jin 



(29) 



whereas k denotes the order of the iterative steps. 

In the case of the memristor bridge circuit, x(t) is the 
normed extent of the space-charge region. The iteration should 
be performed without considering the window function. The 
window function suppresses changes of the state equation 
within the marginal area of the memristor. Therefore for xi' +1 ^ 
the argument of the window function would be of the same 
iteration step. This would lead to a recursive mathematical 
expression and shows one weakness of the Picard Iteration. 
The iteration should be only performed for higher frequencies, 
thus, the window function can be neglected. By setting the 
initial conditions it should be noted that the internal states 
stay within the domain. Attention should be paid to the case 
x > 1 because the memristance would become negative. For 
memristive systems using the HP memristor model, the Picard 
iteration could be performed by 



x\ + ' '-—xq- 



, D 2 



v W 
M(x^) 



■dt , te {t ,t + e} (30) 



with V M lf ] = V Un (t,xf ] 



xjj, ...,Xn S ). Similar to chapter |V 



the used notation for the memristance M(x„ ) is M}^ with 
n,teft. NumW and Den'*' should established as well. The 
only difference to their analogons in chapter [V] is the usage 
of M(xn) instead of M(x„). 

1st iteration step: Using the initial states x„ =x n (0), = 
as arguments for the first iteration step, then the equations for 
the currents are to be 



w/ 1 =v s . 



Im J =V s - 



Im| 1] =V S - 



I M i J =V S - 



M 



1 / NumW+M^-M^-RjA 

1 Kt-^ (31) 



DenPl 



1 NumW+MJ, 11 -Ml. 11 -R l 



Mk 



DenW 



! /, NumW+Mp'-M^-RL' 



1- 



Mk 



DenW 



Mr!' 1 



1 NumW+M^-M^-RL 
DenW 



(32) 



(33) 



(34) 



Mr. 



1] 



2nd iteration step: 

d*j? 
At 



±const ■ Vq ■ sin(G) • t) ■ Mr„ 



(35) 
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By dividing of the variables this equation is 



Ax = ±const ■ Vq ■ M R l „ 1 ■ / sin(o) • t )dt (36) 
Jo 

and leads by integration to 

,P] = ^] ± £^. MR [i]. (1 _ cos(a) . 0) (37) 

[2] [21 

whereat for x\ ' 4 the sign is "+" and for x l 2 3 the sign is negative. 
For this circuit for the second step the Picard Iteration gives a 
good quality solution (regarding Fig.|VI[). Therefore no further 
iteration steps should be performed. The voltage over the load 
resistor can be solved by 

Rl (M®-M®-M? ] -M®). (38) 



v RI J 2 ] 



DenPl 




numerical 

analytica l d as hed line: numerical 

solution. 

Red solid line: analytical 
solution by Picard Iteration 
(second iteration step). 
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Fig. 8. Vrl over time t, analytical and numerical solution. 



VII. CONCLUSIONS 

At this paper a circuit consisting of HP memristors with 
frequency selective behavior is presented. In contrast to higher 
frequencies for low frequencies there is a rectifier like behav- 
ior. This functional change with frequency is reasoned by the 
delay of changing of the internal state caused by an external 
source. To estimate the dynamical behavior in circuits, the 
time behavior of a single memristor in series with a periodic 
source was investigated. The configuration of the circuit leads 
to a nonlinear system of differential equations which describes 
the internal states. Using the Picard Iteration is one possibility 
to solve this system analytically. The frequency selective be- 
havior can be used to realize two modes in one circuit (Graetz 
and Wheatstone circuit). Another applications are a frequency 
controlled regulator or a programmable synaptic membrane 
voltage generator for Spike-Time-Dependent-Plasticity. 
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